Position and frequency shifts induced by massive modes of the gravitational wave 

background in alternative gravity 
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(GR). The deflection and frequency shift effects on light from distant sources propagating through a 
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PACS numbers: 04.30, 04.30.Nk, 04.50. +h, 98.70.Vc 



I. INTRODUCTION 

Einstein's theory of GR has been tested in its weak- 
field approximation and found to pass all the available 
experiments at terrestrial and Solar System scales \vL 
Outside the Solar System, the binary pulsar 1913+16 [2] 
provides indirect evidence for gravitational waves with an 
energy loss consistent with GR. However, other theories 
may produce the same change of orbital parameters: for 
example, the emission of scalar radiation from this binary 
within the context of scalar-tensor gravity is necessarily 
small because, due to the high symmetry of this system, 
the dipole moment is small. As a consequence, the con- 
straints on scalar-tensor gravity imposed by the binary 
pulsar are not competitive with those from Solar System 
tests (however, the binary pulsar data are sufficient to 
rule out Rosen's bimetric theory) []]]. 

Gravitational lensing has provided evidence for light 
deflection on galactic and cluster scales but, due to our 
ignorance of the detailed mass distribution of the lens, 
gravitational lens systems constitute poor tests for the 
theory of gravity (even assuming the validity of GR, the 
lens model is not unique). Instead, one tries to obtain 
information about the mass distribution in the lens by 
assuming the validity of GR and, in this context, obtains 
evidence for dark matter. 

No deviations from Einstein's gravity have been de- 
tected so far in the Solar System or binary pulsar and 
therefore, from the experimental point of view, there is 
no compelling reason to study alternative gravity theo- 
ries. On the other hand, high energy theories that incor- 
porate gravity, such as superstring theory, supergravity, 
and braneworld models, predict deviations from GR in 
the form of extra scalar, vector, or tensor fields of gravita- 
tional origin, massive gravitons, large extra dimensions, 
higher order corrections to the Einstein equations, or vi- 
olations of the Equivalence Principle. The low-energy 



limit of these theories resembles more scalar-tensor or 
f(R) gravity than GR 0, 0]- This fact, in itself, con- 
stitutes a motivation to explore astrophysical and other 
effects in gravitational theories beyond GR. Further, the 
1998 discovery [H] that, if GR is correct, 75% of the en- 
ergy content of the universe is in a mysterious and exotic 
form called dark energy, which propels the accelerated 
expansion of the universe in the present era 0|, leads one 
to be more inclined towards exploring alternative theo- 
ries of gravity rather than reinforcing one's faith in Ein- 
stein's theory which, after all, has been tested only at the 
post-Newtonian level and mostly at Solar System scales. 
It is true that the backreaction of local inhomogeneities 
in an otherwise Friedmann-Lemaitre-Robertson- Walker 
universe certainly affects its dynamics, and that this ef- 
fect is obtained in pure GR without advocating dark en- 
ergy or modified gravity 0. However interesting this 
possibility may be, it has not been possible to produce 
evidence that the magnitude of this backreaction effect is 
such that it can explain the cosmic acceleration observed. 
Backreaction, dark energy, and modified gravity are still 
open possibilities, each scenario has its own difficulties 
and, at present no choice between them is compelling, 
and can be motivated other than by aesthetical consider- 
ations or taste. Of course, a conservative relativist could 
argue that since no deviation from GR has ever been de- 
tected, it is pointless to actively pursue competing theo- 
ries of gravity. On the other hand, proponents of the high 
energy physics point of view would be justified in replying 
that we may actually be detecting the first (large scale) 
deviations from Einstein's theory in the cosmic accelera- 
tion, and that it would be foolish to ignore them. 

We remind the reader that the need to postulate dark 
matter in order to explain the rotation curves at galac- 
tic and cluster scales, has led people to doubt not only 
general relativity, but even Newtonian gravity, and has 
produced MOND and TeVeS theories [& EJ, and f(R) 
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gravity with anomalous couplings to matter [lfjj]. There- 
fore, it seems reasonable to try to identify possible ways 
to test gravity beyond the Solar System. There seem to 
be two conflicting points of view, corresponding to two 
different communities: deviations from Einstein's gravity 
are regarded as unavoidable by high energy physicists, 
their detection being only a matter of technological lim- 
its. On the other hand, classical relativists may regard 
the theories producing such deviations as exotica, and it 
is certainly true that the latter do not have experimental 
support so far. While theories of gravity alternative to 
general relativity are purely hypothetical so far, they are 
theoretically well-motivated and there is scope to try to 
detect deviations from Einstein's gravity. It is interesting 
to point out that preliminary results in positive energy 
theorems exist for such theories as discussed in details in 
Besides, in the light of developments in high energy 
physics, such possibilities should not be discarded a pri- 
ori. Eventually, experiment is the judge and the failure 
to detect deviations from GR further constrains alterna- 
tive theories and worsens the fine-tuning problems that 
they may have. 

Here we do not want to argue in favour of GR or its 
competitors: rather, we try to bridge the two points of 
view and we study possible deviations from GR predicted 
by high energy theories. We focus on a possible astro- 
physical effect that was studied in the past in the context 
of GR, and found to be negligible, but is potentially inter- 
esting in alternative theories of gravity. This effect con- 
sists of the deflection and frequency shift of a light beam 
due to its propagation through a stochastic background 
of gravitational waves. In many theories of gravity, ex- 
tra gravitational fields (scalar, vector, and tensor) appear 
in addition to the usual massless spin 2 graviton famil- 
iar from GR. These modes, massless or massive, corre- 
spond to extra degrees of freedom contained in the metric 
tensor and show up as gravitational waves emitted 
by early astrophysical sources or excited by cosmological 
processes, superposing to form a stochastic background. 
Such a background, analogous to the cosmic microwave 
background of electromagnetic waves, is well known in 
GR, and the propagation of light rays through it has been 
studied in detail El EE El EE EE 

Il7| . Consider a pencil 
of light rays propagating from a distant light source (pos- 
sibly at a cosmological distance) to an observer. Since 
gravitational waves deflect light rays and perturb their 
frequency, naively one expects a photon undergoing N 
scatterings in this background to be described by a ran- 
dom walk and its deflections, or frequency shifts, to add 
stochastically as y/N. Even though the deflection is at 
most of linear order in the gravitational wave amplitudes 
which are very small, since the travelled distance can be 
large, such a cumulative (or "L-") effect that grows as \J~L 
could compensate for it, and it has indeed been claimed 
in the past EE EH- Intuition fails, however, because it is 
based on familiarity with random walk processes in which 



the scatterers are static or nearly static, while the mass- 
less gravitons of GR, responsible for photon scattering, 
propagate at the speed of light. The size of the deflection 
(or frequency shift) effect is a matter of relative velocities, 
i.e., of the difference between the speed of the propagat- 
ing signal and that of the perturbations from a uniform 
background through which the signal propagates. When 
this fact is taken into account, the cumulative L-effect 
disappears 3, IE IE 12, Ji|- The quantitative descrip- 
tion of the deflection (or frequency shift) effect depends 
not only on the relative speed, but also on the spin s 
of the field responsible for the non-stationary perturba- 
tions in the otherwise homogeneous medium. A com- 
prehensive quantitative treatment is given in [lfij |. The 
analogous situation for massless scalar modes in scalar- 
tensor gravity was briefly considered in (2(| and it was 
found that, in spite of a logarithmic dependence of the 
rms deflection on L, the effect is numerically comparable 
to the one in GR and, therefore, completely negligible 
for practical purposes. However, the spectrum of gravi- 
tational theories now available is considerably larger and 
the consideration of astrophysical effects due to massive 
fields of various spins forming a stochastic background 
can potentially be of interest, since massive fields can al- 
low for a cumulative L-effect, which will be explored in 
the following sections. 

The plan of this paper is as follows. In Sec. II we 
briefly recall the physics of the deflection and frequency 
shift effects for non-stationary perturbations of different 
spins. In Sec. Ill this analysis is applied to gravitational 
theories that predict deviations from GR. In Sec. IV the 
case of modified (or f(R)) gravity is studied in detail, 
while Sec. V contains a discussion and the conclusions. 



II. DEFLECTIONS AND FREQUENCY SHIFTS 
CAUSED BY PROPAGATION IN A 
GRAVITATIONAL WAVE BACKGROUND 



To realize how gravitational waves induce deflections 
and frequency shifts in a light ray with tangent p^ 
that traverses them, it is sufficient to consider the null 
geodesic equation 



(1) 



By locally expanding the metric as <? M „ = rj^ + h^ v in an 
asymptotically Cartesian coordinate system, where the 
perturbations h^u (with |/i, M „| << 1) describe gravita- 
tional waves, computing the Christoffel symbols Y^g to 



first order, and using the fact that p M = p'Ls 



Spv 



(1,0, 0, 1) + Sp^ with <5p M = 0(h) for a photon with un- 
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perturbed path along the z-axis, one obtains 

i-O 



d ^ r a/3 P(0)P(0) 



dz (/ioo ~ 2/i 03 + W + 0(h 2 ) , (2) 



While Linder, in the context of GR, focussed on mass- 
less spin 2 gravitons and the limit v — > 1, here we are 
interested in the opposite limit for massive modes. As 
shown in the next section, some of these modes can be- 
come very massive, corresponding to v — > 0. In this case 
J_i becomes negligible and we are left with the J_2 con- 
tribution. 



where the integral is computed along the unperturbed 
path from the source S to the observer O. This 
shows that, in GR, a gravitational wave propagating 
(anti) parallel to the light ray has no effect on it, to first 
order. 1 If the describe a superposition of many 
waves with random phases, directions of propagation, 
and polarizations, one will obtain deflections such that 
(<5p M ) = but ((Sp^) 2 ) ^ 0. Therefore the problem is 
whether these random deflections (for /i = 1,2,3) or fre- 
quency shifts (for fi = 0) add stochastically. This prob- 
lem has been solved by Linder in a more general 
context by considering random fluctuations due to inho- 
mogeneities propagating with arbitrary speed v between 
the light source (at z — 0) and an observer (at z = L) 
and due to a superposition of fields of spin s = 0, 1, or 2. 
By writing the deflection due to a single mode as 



dz e 



(3) 



and e(t,x) = TZe (eoe 1 *^ 2 ^) , Linder obtains the mean 
square deflection 



-is 2 



2 s 



=0 (Tle 2 e s ) ^ a nJ n 



(4) 



n=-2 



where s is the spin of the field responsible for the inho- 
mogeneities, a n are constants, and J n are the integrals 



Jn 



1 



(kL) 



n+l 



dy y n sin 2 y 



(5) 



One is interested in the limit for wavenumbers k and 
lenghts L such that kL >> 1; in this limit the integrals 
J n for n > cannot cause an L-effect and we focus on 
the integrals for n = — 2, — 1, given by 



J_i = -4 (1 + s) « (1 - v 2 



2 , sin y , „ 

dy , 6 

kL(.i+v) y 



J-2 = kL[l-v 



2\!+ s 



, sin y 

dy ~~2 ■ (7) 



III. APPLICATION TO ALTERNATIVE 
THEORIES OF GRAVITY 

In several alternative theories of gravity, massive grav- 
itational fields appear which can potentially give rise to 
an L-effect. Some of them are inspected in the following. 



III.l Scenarios with large extra dimensions 



It has been suggested [2l| that the hierarchy prob- 
lem could be solved in theories with large (sub-millimeter 
size) extra spatial dimensions, in which gravitons prop- 
agate through (3 + n)-dimensional space while non- 
gravitational physics is confined to the ordinary three 
spatial dimensions (see [22] for a review). The n ex- 
tra dimensions are compactified, e.g., on a torus with a 
radius R n and gravity can be strong already at the TeV 
scale. The gravitons propagating in the extra dimensions 
acquire a mass given by 



2 47T71 2 



where 



32-17n 

Rn = 2 • 10^^ cm 



(«) 



(9) 



The model is ruled out for n = 1 and marginally ruled 
out for n — 2 (for which Ri ~ 2 mm), but is viable 
for n > 2, corresponding to R n < 10~ 6 cm. The dis- 
persion relation fc M fc M = — m 2 for the massive gravi- 



tons yields the group velocity v g — 



If 



9 ~ y/ml+k 2 

is sufficiently large, many (most) modes composing the 
gravitational wave background will have k « m n and 
v g ~ ck/rrin << c. For example, for n — 3, one ob- 
tains m3 ~ 3 • 10~ 19 / 3 cm' 1 ; for waves of wavelength 
X g ~ 10 3 km it is v g ~ 10~ 2 c, while longer waves with 
A 9 

and X g ~ 10 3 km, it is v g ~ 10 



3 • 10 s km= 2 A.U. yield v g - 10~ 7 c. For 77 = 4 



-17 



C. 



III.2 / (R, R^R^ , R^^R^i" 7 ,n\R,a 2 R, ...) theories 



1 This can be seen by adopting the transverse-traceless gauge in 
which /ioo = ^03 = ^33 = for a gravitational wave propagating 
in the ±2 direction. 



In general, in theories described 

by a Lagrangian density of the type 
f (R,R^ u R' J - v ,R pi/ p (r R' MV P' T ,nR,n 2 R, ...), there are 
scalar, vector, and tensor modes, massive or massless, 
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and these can, in principle, contribute to the grav- 
itational wave background and produce an L-effect. 
However, some of these massive modes are ghosts, which 
precludes further consideration of these theories. An 
exception are theories with Lagrangian of the form 
/ (R, G), where Q = R 2 - AR^W"' + R^ pa R^ pa is the 
Gauss-Bonnet combination. At least if certain conditions 
are satisfied, ghosts are avoided in these theories [Hj]. 



III. 3 N = 2, 8 extended supergravity 

The supergravity multiplet in N — 2, 8 extended su- 
pergravity contains a graviton, a gravivector field, two 
Majorana gravitinos for N = 2, and a graviscalar field for 
A = 8. The graviscalar violates the Weak Equivalence 
Principle j24j], and both graviscalar and gravivector are 
short-ranged. The available experiments set the limits 
on their ranges Ri and R a , respectively, 24, IB] 



Ri < 0.6 cm, Ri > 13 cm (N = 2), (10) 



Ri < 0.4 cm, Ri > 40 m (N = 



(11) 



R a < 0.15 cm, 60 m < R a < 100 m. (12) 

If these fields are truly short-ranged, they can also con- 
tribute as massive modes to the gravitational wave back- 
ground and the analysis of the previous section applies. 

There is scope, therefore, to consider the limit v — ► 
for massive gravitons in these scenarios in the discussion 
of the previous section. 



III. 4 rms deflections and frequency shifts due to 
massive modes 



Since J_i — » in the limit of heavy modes v 
are left with the contribution of 



kL 



sin 2 y 



+kL/2 

dy , 

-kL/2 V 



in eq. (0]). By using 



dy 



sin 2 y 



cos(2y) 2ySi(2y) - 1 



2y 



2y 



0, we 



(13) 



(14) 



where Si(z) = dt ^ = | - J+°° dt ^ is the sine 
integral, one obtains 



J_ 2 = 2 [cos(fcL) + kLSi{kL) - 1] 



(15) 



in the v — > limit. The term 2 [cos(fcL) — 1] assumes val- 
ues in the interval [—4, 0] and oscillates as kL becomes 
large, while the second term kLSi(kL) dominates. Since 



Si(+oc) 
tion 



7r/2, the limit kL >> 1 yields the rms deflec- 



\ 



a- 2 



■nkL 



(16) 



s=0 



for these modes, where an L-effect is indeed present and 
can, in principle, compensate for small values of the grav- 
itational wave amplitudes £2 to produce a non-negligible 
effect. This is not surprising since in the limit v — > the 
propagation of the photon reduces to a random walk. 

More precisely, keeping the dependence of J_2 on v 
yields 



J_ 2 ~ 2 (1 - v 2 ) kL [Si{kL) + vSi(kLv)} 



(17) 



and 



4k) 



. ^2a- 2 kL{\ 

\ s=0 



v 2 )[Si(kL)+vSi(kLv)} (Tie 2 



(18) 



In all the scenarios listed above one can expect very 
massive modes for which an L-effect exists and the rms 
deflection or frequency shift is given, in order of magni- 
tude, by 



'kLe 



(19) 



where e is the magnitude of the wave amplitude for the 
massive mode considered. The estimation of this quan- 
tity is difficult because it depends on the processes gen- 
erating the cosmological background, which are subject 
to much speculation and large uncertainties even in GR. 
The calculation of precise spectra of gravitational modes 
in specific processes is beyond the purpose of this work. 
We assume that detailed studies can provide, in princi- 
ple, estimates of e in various frequency bands following 
assumptions about specific generating processes. In this 
paper we study in more detail the case of f(R) gravity. 



IV. MASSIVE MODES IN f(R) GRAVITY 

Modified or "/(i?)" gravity has been proposed recently 
in order to explain the current acceleration of the uni- 
verse without resorting to dark energy [H, [27I . f(R) 
gravity comes in three versions: the metric [26ll27|. Pala- 
tini [2^, and metric-affine [2^] formalisms. In the metric 
formalism, in which the metric tensor is the only indepen- 
dent variable and the connection is the metric connection, 
the action is 



A = 2l J d ' x ^~9tt R ) + SimaUer) > 



(20) 
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where f(R) is a non-linear function of its argument re- 
placing^ the usual Emstein-Hilbert Lagrangian R — 2 A 
([26, |T 



7j — see [30j, |31|, |32| for reviews). Corrections to 
this Lagrangian that become important as R — > can 
explain the current acceleration of the universe without 
resorting to dark energy, while early universe physics in a 
strong curvature regime is instead affected by corrections 
described by positive powers of R. Indeed, the renormal- 
ization of GR introduces quadratic corrections [33] , a fact 
that was exploited in Starobinsky's scenario of inflation 
without scalar fields [34j]. The condition f"(R) > is re- 
quired in the metric (but not jn the Palatini) formalism 
for the absence of tachyons 
stability 



35, 36] and for non-linear 



Metric f(R) gravity is dynamically equivalent to an 

with a non-trivial po- 
: f'(R), an equivalent 



uj = Brans-Dicke theory [31 
tential. In fact, by setting 
action is 



S=^- f d 4 x^g~ [cj>R - VU)} + S {matter) , (21) 
2k J 

where 

V(cf>) = 4>R(<P) - f(R(4>)) . (22) 
The scalar degree of freedom f'(R) satisfies the equation 



from which one obtains the effective mass 131] 



, Rf"{R)-f'{R) 



3f"(R) 



(23) 



(24) 



with R = R(<p). 

In metric f(R) cosmology, the dependence of the ef- 
fective mass of 4> 011 the curvature and, therefore, on 
the environmental density is exploited in the chameleon 
mechanism in order to make these theories viable. At 
Solar System densities, the scalar has a very short 
range, thus evading the constraints imposed by Solar 
System and terrestrial experiments on the equivalent 
Brans-Dicke theory, while at cosmological densities this 
range becomes very long and can affect cosmology. This 
chameleon mechanism (well-known in quintessence mod- 
els [1^]) makes these theories viable, but at the same 
time it renders the long wavelength scalar modes forming 
the stochastic background effectively massless. There- 
fore, the analysis of massless Brans-Dicke scalar modes 
of Ref. [2(3] applies and no L-effect is present. 

It is more interesting, from this point of view, to con- 
sider f(R) theories relevant for early universe physics. 
For example, in the model f(R) = R + aR 2 , it is 
Tiiffr — 1/V6a. One expects the parameter a weight- 
ing quantum corrections to the Einstein-Hilbert action 



to be small and, hence, a large mass for the scalar de- 
gree of freedom <f>, which propagates with group velocity 
v g ~ cfc/m^ = \/6acfc. 



Assuming the conformal transformation 



2$ 



with 



„2<I> 



f'{R) , (25) 



where the prime indicates differentiation with respect to 
the Ricci scalar R and <f> is the "conformal scalar field", we 
obtain the conformally equivalent Hilbert-Einstein action 



A ~ 2k 



x\j —g 



i? + £ ($,<&;„) 



(26) 



where L (<£>,<I>; At ) is the conformal scalar field contribution 
derived from 



Ruv 



and 



g^ ;S 



R 



-2<S> 



R - 6D$ - 6$.a$ ; ' 



(27) 



(28) 



In any case, as we will see, the C ( < f>,$- AI )-term does not 
affect the gravitational wave tensor equations so it will 
not be considered further. 2 

Beginning with the action (f26|) and deriving the 
Einstein-like conformal equations, the gravitational wave 
equations expressed in the conformal metric are 



□ hi = 



(29) 



Since no scalar perturbation couples to the tensor part 
of the gravitational waves, we have (5$ = and then 



3 



e- 2 V'e 2 * 



Sgu = h{ 



(30) 



which means that h\ is a conformal invariant. As a con- 
sequence, the plane- wave amplitude defined by h\ (t, x) — 
h(t) ej exp(ifcix'), where e\ is the polarization tensor, are 
the same in both metrics. In any case, the d'Alembert 
operator transforms as 



□ = 



-2* 



(□ + 2$ ;A V; 



(31) 



and this means that the background is changing while 
the tensor wave amplitude is not. 

In order to study the cosmological stochastic back- 
ground, the operator lf3lj) can be specified for a 
Friedmann-Robertson- Walker (FRW) metric given by 



-dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) 



(32) 



2 Actually, a scalar component of gravitational radiation is often 
considered [13, S3], but here we are taking into account only the 
genuine tensor part of the stochastic background. 
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and then eq. (|29| becomes 

h + (3H + 2$) h + k 2 a- 2 h = , (33) 
d 2 d 

where □ = — — + 3H— and k is the wave number. 
ot z at 

It is worth stressing that eq. ([33)1 applies to any f(R) 
theory whose conformal transformation can be defined 
as e 2 * = f'(R). The solution, i.e., the gravitational wave 
amplitude, depends on the specific cosmological back- 
ground (i.e. l _a(t)) and the specific theory of gravity 
(i.e., $(t)) [42|. Considering also the conformal time 
drj = dt/a, eq. lj33"|) reads 



d 2 h 2 d\ dh 
drj 2 x dr/ drj 







(34) 



where x = oe*. Inflation means that a(t) — a exp(Ht) 
and then T) = J dt/a = (aH)' 1 and -^f- = -r/' 1 . The 



exact solution of (1341) is 



h(rj) = V2 k 2 [Ci sin fcry + C 2 cos fcij] 



(35) 



Inside the H^ 1 radius we have krj ^> 1. Furthermore, 
considering the absence of gravitons in the initial vacuum 
state, we have only negative-frequency modes and then 
the adiabatic behavior is 



2 fc 1/2 — Cexp(-ikrj) . 
7r aH 



(36) 



At the first horizon crossing (aH — k), the averaged 
amplitude of the perturbation Ah — (k/2ir) 3 / 2 \h\ is 



A h 



C 
2^2 



(37) 



When the scale a/k grows larger than the Hubble radius 
H~ x , the growing mode of evolution is frozen, that is, 
it is constant. This situation corresponds to the limit 
krj -C 1 in eq. (|35| . Since $ acts as the inflaton field, 
it is $ ~ at re-entry after the end of inflation. Then 
the amplitude Ah of the wave is preserved until the sec- 
ond horizon crossing after which it can be observed, in 
principle, as an anisotropy perturbation in the cosmic mi- 
crowave background. It can be shown that AT/T < Ah 
as an upper limit to Ah since other effects can contribute 
to the background anisotropy [13]. From these consid- 
erations, it is clear that the only relevant quantity is 
the initial amplitude C in eq. lf36"ll . which is conserved 
until re-entry into the horizon. Such an amplitude di- 
rectly depends on the fundamental mechanism generat- 
ing the perturbations. Inflation gives rise to processes 
capable of producing perturbations as zero-point energy 
fluctuations. Such a mechanism depends on the theory 
of gravitation adopted and then (AT/T) could consti- 
tute a further constraint to select a suitable /(-R)-theory. 



Considering a single graviton in the form of a monochro- 
matic wave, its zero-point amplitude is derived through 
the equal time commutation relations 



[h(t, x), Tr h (t, y)] = i 5 3 (x - y) 



(38) 



where the amplitude h is the field and Hh is the conjugate 
momentum operator. Writing the Lagrangian for h 



(39) 



in the conformal FRW metric g^ v (h is conformally in- 
variant), we obtain 



■Kb, = 



dC 
dh 



e 2 *a 3 A 



Then, eq. <[38|l becomes 

h(t, x), h(t, y) 



6 3 (x~y) 



(40) 



(41) 



and the fields h and h can be expanded in terms of cre- 
ation and annihilation operators 



h(t, x) 



h(t, x) 



1 



(2tt) 3 / 2 



1 



d 3 k \h(t)e 



-ikx 



h*(t)e 



-ikxl 



d A k 



h(t)e 



—ikx 



h*(t)e 



-\-ikx 



(42) 



(43) 



(2tt) 3 / 2 

The commutation relations in conformal time are then 

i(2n) 3 



[hti* - h*h'} = 



(44) 



The substitution of eqs. (J36J and (J3ZJ) yields C = 
V2Tv 2 He~*, where H and $ are calculated at the first 
horizon-crossing and then 



(45) 



which means that the amplitude of gravitational waves 
produced during inflation directly depends on the given 
f(R) theory since $ = \ In f'(R). Explicitly, it is (42 1 



Ah 



H 



(46) 



where f'(R) > is necessary in order for the graviton 
to carry positive kinetic energy [3l|. The representation 
of f(R) gravity as a Brans-Dicke theory is particularly 
useful when dealing with the scalar component of gravi- 
tational waves, ruled by the equation [44| 



□$ = TO $ 



(47) 



where $ = —S4>/<po- The scalar field generates a third 
component for the tensor polarization of gravitational 
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waves and the total perturbation describing a gravita- 
tional wave propagating in the positive z direction is 



the redshift at the equivalence epoch has to be consid- 
ered. Taking into account Ref. [49], one gets 



h^(t -z) = A + (t - z) e ( +J +A x (t-z)e 



(x) 



Qsgw (/*) — 



16 



PdS 



9 ppi 



ancfc 



-(i + z eq y 



(55) 



-$(t- 



(48) 



The term A+(t- z)e$ +A X (t 
standard (i.e., tensorial) polarizations of a gravitational 
wave arising from GR in the TT gauge [45j], while the 
term $(i — z)efy is the extension of the TT gauge mode 
to the scalar case. Three different degrees of freedom are 
present (see eq.(32) of jilj), while only two are present 
in standard GR. Then, for a purely scalar gravitational 
wave, the metric perturbation is [4J| 



- z e 



(x) 

/UJ 



for the waves which, at the epoch in which the universe 
becomes matter-dominated, have a frequency higher than 
describes the two H eq , the Hubble parameter at equivalence. This situa- 



(49) 



The stochastic background of scalar gravitational waves 
can be described in terms of the scalar field $ and char- 
acterized by a dimensionless spectrum (see the analogous 



definition for tensor modes in 140, 46, 47, 4 



^sgw (/*) 



1 dp 



where 



Pc d\nf 



ZH 2 



8ttG 



(50) 



(51) 



is the (present) critical energy density of the universe, 
Hq is the Hubble parameter today, and dp sgw is the en- 
ergy density of the scalar gravitational radiation in the 
frequency interval (f,f + df). We are now using stan- 
dard units. Now it is possible to write an expression for 
the energy density of the stochastic scalar relic gravitons 
background in the angular frequency interval (uj, uj + du) 
as 



dp, 



sgw 



2hu> 



( UJ 2 dbJ 



V 2tt 2 c 3 



N„ = 



HH 2 dS H 2 duj _ HH 2 dS H 2 df ^ 

47T 2 C 3 LJ 4lT 2 C 3 / ' 

(52) 

where /, as above, is the frequency in standard comoving 
time. Eq. lf5"2"j) can be rewritten in terms of the critical 
and de Sitter energy densities 



m = 



8irG Pc 



Has = 



inGpds 
3c 2 



Introducing the Planck density ppi a 
trum is given by 



KG 2 



(53) 



, the spec- 



^sgw (/*) 



1 dp a 



f dp s 



Pcdlnf p c df 



16 PdS 

9 ppianck 



(54) 



At this point, some comments are in order. First, the 
calculation works for a simplified model that does not in- 
clude the matter-dominated era. If the latter is included, 



'e<j, 

tion corresponds to frequencies / > (1 + z eq yl 2 H.Q to- 
day. The redshift correction in eq. |55|) is needed since 
the present value of the Hubble parameter H would be 
different without a matter-dominated contribution. At 
lower frequencies, the spectrum is given by 4^, 4^| 



ttsgw(f) oc / 



(56) 



As a further consideration, let us note that the results 
(|54|l and (|55|l . which are frequency-independent, do not 
hold in the entire range of physical frequencies. For waves 
with frequencies less than the present Hubble parameter 
Ho, the notion of energy density is not defined because 
the wavelength becomes longer than the Hubble scale. 
Similarly, at high frequencies, there is a maximal fre- 
quency above which the spectrum rapidly drops to zero. 
In the above calculation, the simplifying assumption that 
the phase transition from the inflationary to the radia- 
tion dominated epoch is instantaneous has been made. 
In the physical universe, this process occurs over some 
time scale At, with 



fmax — 



a(h) 1 
a(t a ) At 



(57) 



which is the redshifted rate of the transition. In any 
case, flsgw drops rapidly. The two cutoffs at low and 
high frequencies for the spectrum guarantee that the total 
energy density of the relic scalar gravitons is finite. For 
GUT-scale inflation, it is of the order [46] 



Pds 



PPlanck 



1(T 



(58) 



These results can be quantitatively constrained consider- 
ing the recent WMAP release. In fact, it is well known 
that WMAP observations put severe restrictions on the 
spectrum. In Fig. [TJ the spectrum fl sgw is mapped: con- 
sidering the ratio pds/ ppianck, the relic scalar gravita- 
tional wave spectrum seems consistent with the WMAP 
constraints on scalar perturbations. Nevertheless, since 
the spectrum falls off as f~ 2 at low frequencies, today 
at LIGO '/VIRGO and LISA frequencies (indicated in 
Fig. [TJ, one gets 



n sgw (f)h 2 oa < 2.3 x 10" 



(59) 



where hwo = Ho/ (100 km • s _1 • Mpc -1 ). It is interest- 
ing to calculate the corresponding strain at / ~ 100Hz, 
where interferometers such as VIRGO and LIGO achieve 
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maximum sensitivity. The well known equation for the 
characteristic amplitude [iiL |48| adapted to the scalar 
component of gravitational waves 



$ c (/)~1.26x 10 



-18 



1H2 

~T 



v^yu/i , (60) 



can be used to obtain 

$ c (100Hz) < 2 ■ 10~ 26 



(61) 



Then, since we expect a sensitivity of the order of 10~ 22 
for the above interferometers at / ~ 100 Hz, we need 
to gain four orders of magnitude. Let us analyze the 
situation also at lower frequencies. The sensitivity of 
the VIRGO interferometer is of the order of 10" 
/ ~ 10 Hz and in that case it is 

$ c (10 Hz) < 2 • 10~ 25 . 



at 



(62) 



The sensitivity of the LISA interferometer will be of the 
order of 10 -22 at / ~ 10~ 3 Hz and in this case it is 



$ r (10- J Hz) < 2 ■ 10- 



(63) 



This means that a stochastic background of relic scalar 
gravitational waves could, in principle, be detected by 
the LISA interferometer. 



LISA LIGO 



15 -10 -5 



Hz {- Graphics -> 



Figure 1: The spectrum of relic scalar gravitational waves in 
inflationary models is flat over a wide range of frequencies. 
The horizontal axis is log 10 of frequency, in Hz. The vertical 
axis is log 10 D. as -w The inflationary spectrum rises quickly at 
low frequencies (wave which re-entered the Hubble sphere af- 
ter the universe became matter-dominated) and falls off above 
the (appropriately redshifted) frequency scale f max associ- 
ated with the fastest characteristic time of the phase transi- 
tion at the end of inflation. The amplitude of the flat region 
depends only on the energy density during the inflationary 
stage; we have chosen the largest amplitude consistent with 
the WMAP constraints on scalar perturbations. This means 
that, at LIGO and LISA frequencies, Q 3 gw < 2.3 ■ 1CT 12 . 



To estimate the rms deflection effect from massive 
scalar modes in the gravitational wave background, we 
restrict to periods that are of the order of hours or days. 



A longer period would result in a "frozen" effect which 
is much less likely to be detected, while a much shorter 
period would probably render the deflections unobserv- 
able because only an averaged position shift would be 
recorded during observation times longer than the pe- 
riod itself and a slightly blurred image would be the out- 
come (although fast photometry might allow to push the 
limits). Assuming a frequency ~ 10 -3 Hz at a distance 
L ~ 500 kpc and using the upper limit (|63"1) . one obtains 
a rms deflection 9 rms ~ VkLe ~ 10~ 10 . The maximum 
resolution expected with high precision astrometry is of 
the order of microarcseconds (~ 10 -7 radians), three or- 
ders of magnitude above the required sensitivity for de- 
tection. For galactic sources at L ~ 5 kpc, to which high 
precision astrometry is more likely to apply, rms drops 
by another order of magnitude. One can, of course, con- 
sider different sources of electromagnetic radiation with 
slightly higher frequency, at more promising distances L, 
and perhaps find mechanism which produce higher scalar 
amplitudes <&: at a first look, however, it is unlikely that 
the four orders of magnitude necessary for detection can 
be bridged in the foreseeable future. 



V. OUTLOOKS 

All modern theories of high energy physics unifying 
gravity with the other interactions predict departures 
from GR; however, no such deviation has been observed 
so far in Solar System experiments, and practically all the 
experimental constraints on such deviations are obtained 
within the Solar System (the binary pulsar and grav- 
itational lensing provide constraints that are not com- 
petitive with those obtained from Solar System experi- 
ments). Therefore, it is interesting to explore astrophys- 
ical effects outside of this narrow region of the universe 
that could potentially exhibit deviations from Einstein's 
theory. The cumulative deflections, or frequency shifts 
due to propagation of light from distant sources through 
random massive modes of the gravitational wave back- 
ground could constitute such an effect. There is now a 
wide range of theories predicting massive scalar, vector, 
and tensor modes that can lead to such an effect. How- 
ever, the astrophysical and cosmological processes gen- 
erating cosmological gravitational wave backgrounds in 
these theories are still unexplored. Here we do not at- 
tempt to estimate the average strength e of the various 
modes appearing in these theories of gravity, in different 
ranges of wavelengths, and under various assumptions. 
We limit ourselves to outline an estimate for a particular 
case: massive scalar modes in f(R) gravity. This class 
of theories has been the subject of much recent atten- 
tion in order to explain the observed acceleration of the 
universe without resorting to dark energy. For the situ- 
ation considered here, the deflections seem to be a few 
orders of magnitude too small for detection; on the other 
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hand, we find that these scalar modes are certainly of 
interest for direct attempts at detection with the LISA 
experiment. According to our preliminary discussion, the 
indirect detection with the position and frequency shift 
effect does not seem to be feasible with current technol- 
ogy; however, a more detailed analysis is necessary before 
definitive conclusions can be drawn. It is, in principle, 
possible that massive gravitational wave modes could be 
produced in more significant quantities in cosmological 
or early astrophysical processes in alternative theories of 
gravity — the latter are still unexplored. This possibility 
should be kept in mind when looking for a signature dis- 
tinguishing these theories from GR, and seems to deserve 
further investigation. 
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